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ABSTRACT 



The total mass of distant star clusters is often derived from the virial theorem, using line-of-sight velocity dispersion measurements 
and half-light radii. Although most stars form in binary systems, this is mostly ignored when interpreting the observations. The 
components of binary stars exhibit orbital motion, which may increase the measured velocity dispersion, and may therefore result 
in a dynamical mass overestimation. In this paper we quantify the effect of neglecting the binary population on the derivation of the 
dynamical mass of a star cluster. We simulate star clusters numerically, and study the dependence of the derived dynamical mass on 
the properties of the binary population. We find that the presence of binaries plays a crucial role for very sparse clusters with a stellar 
density comparable to that of the field star population (~ 0.1 stars pc -3 ), as the velocity dispersion is fully dominated by the binary 
orbital motion. For such clusters, the dynamical mass may overestimate the true mass by up to an order of magnitude. For very dense 
clusters (3; 10 7 stars pc~ 3 ), binaries do not affect the dynamical mass estimation significantly. For clusters of intermediate density 
(0.1 - 10 7 stars pc~ 3 ), the dynamical mass can be overestimated by 10 - 100%, depending on the properties of the binary population. 

Key words, star clusters: general — methods: numerical — binaries: general 



1. Introduction 

Observations have shown that the majority of the field stars are 
part o f a binary or multiple system (e.g. jDuquennov & Mavor 
Moreover, both observations and numerical simulations 
have indicated that this property is primordial: the vast ma- 
jority of stars a re formed in binar y or multiple systems (e.g. , 
Mathieul fl99H iMason et all 119981: iGoodwin & Kroupal [2 005 ; 



Kobulnicky & Fryerfl2007tlKouwenhoven et al JI2005L I2007L and 
references therein). Although this has been known for more than 
a decade, binaries are often not properly taken into account when 
analysing integrated spectral line data of young star clusters. 

Young star clusters, with typical masses of M c \ = 10 3-6 M Q 
indicate recent or ongoing violent star formation. Their forma- 
tion is often triggered by mergers and close encounters between 
galaxies. Only a fraction of these young massive star clusters 
evolve into old globular clusters, while a substantial fraction 
(~ 60 - 90%) may dissol ve into the field star populat ion within 
about 30 Myr (see, e.g., Ide Griis & Parmentieil 120071, for a re- 
view). In order to study the formation and fate of these star clus- 
ters, it is necessary to obtain good estimates of their total mass, 
stellar content, dynamics, and binary population. In this paper 
we focus on the derivation of the total mass and the properties of 
the binary population in particular. 

There are two straightforward methods to determine the to- 
tal mass of a star cluster. The first is based on a derivation from 
the luminosity of the cluster. This photometric mass determi- 
nation (e.g.. lMengel et alj|2005l: iMa et alj|2006l) is independent 
of the assumptions about the properties of the binary popula- 



tion, as the total mass is derived from the integrated luminos- 
ity, to which each star (whether single, primary, or companion 
star) contributes. This method requires an a priori knowledge 
of the stellar mass distribution /m(M) in the cluster (and hence 
of the mass-to-light ratio), accurate estimates of the age, dis- 
tance, metallicity and interstellar extinction. The second method 
is based on the virial theorem: the dynamical mass, Md yn , is 
derived from the (projected) half-light radius R^i and the line- 
of-sight velocity dispersion, <ti os . The half-light radius is often 
assumed to be equal to the half-mass radius Rh m , i. e., no mass 
segregation is assumed to be present (see, however, Boil v et ail 
2005: lFleck et alj|2006l) . An estimate of M dvn can be obtained 
using the equation derived bv lSpitzeri d 1987b : 



M dyn = rj 



(1) 
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where G is the gravitational constant, and tj is a dimensionless 
proportionality constant. Spitzer's equation is valid under the 
following assumptions: (i) the cluster follows a Plummer den- 
sity model, (ii) all stars are equal-mass stars, (iii) no binary or 
multiple stars are present and (iv) the cluster is in virial equilib- 
rium. An estimate of the accuracy of the derived masses can be 
obtained by comparing photometric and dynamica l masses (e.g., 
i Mandushev et alj Il99lb ISmith & G allagher ] 120011: lLarsen et al ' 
2004t iMaraston et alJl2004t Ide Griis et all l2005f iBastian et al 

120061: lLarsen et al.ll2007l) . 

The Plummer model (PI urnmerj[T91ll) assumed in Eq. (Q]) 
describes the structure and dynamics of mature star clusters 
with reasonable accuracy. It has been in use for a long time 
because of its mathematical simplicity. iKingl d 19621 1 19661) de- 
veloped a set of models, nowadays known as King models, 
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that provi de a more accurate des cription for globular clusters 
(see, e.g.. lMevlan & Heggidll997l) . Although Eq. ([TJ is derived 
for the Plummer model, it is also a good approximation for 
King model s, with a slightly adjusted proportionality constant 
77 (see, e.g.. iFleck et al.l l2006). Very young (^ 1 Myr) clusters, 
however, often exhib it an irregular or flocculent structure (e.g., 
Elmegree n et al.ll2000l) . so that Eq. ([]]) may not be a good ap- 
proximation anymore. The equal-mass assumption is clearly in- 
appropriate, although adopting a more realistic mass distribu- 
tion affects the derived Md yn only mildly (see § I3.2| >. The as- 
sumption of virial equilibrium is probably appropriate for older 
(^ 50 Myr) clusters. It is, however, an incorrect assumption 
for most young (S= 20 Myr) clusters; a su bstantial fraction of 
these may suffer from infant mortal ity (e.g jGoodwin & Ba stian 
2006; Ide Griis & Parmentierll2007l and references therein); we 
will briefly discuss this issue in § 13.11 In this paper, however, we 
focus primarily on the effects of assumption (iii), the assumed 
absence of binary systems, and study how the parameter 77 de- 
pends on the properties of the star cluster and the binary popula- 
tion. 

In a cluster consisting of single stars, the velocity dispersion 
traces the motion of each particle (i.e., star) in the cluster poten- 
tial. In a cluster with binary stars, on the other hand, we do not 
measure the motion of each particle (i.e., the binary centre-of- 
mass), but of the individual binary components. These have an 
additional velocity component due to their orbital motion, which 
may result in an overe s timation of the dynamical cluster mass 
(e.g. JBosch et1ul20inilFle~ck et alJl2006l:lApai et al.ll2007h . For 
a very sparse, dissolving cluster (5 0.1 stars pc~ 3 ), the motion of 
the centre-of-mass of each binary in the cluster is much smaller 
than the orbital motion of the binary components. The spectral 
line width in such a cluster is thus dominated by orbital mo- 
tion, and is not representative of the motion of the binaries in 
the cluster potential. Unless this effect is corrected for, the de- 
rived dynamical mass could be significantly overestimated. On 
the other hand, for a very dense cluster (^ 10 7 stars pc~ 3 ), the ef- 
fect of binaries is almost negligible. Nevertheless, the presence 
of binaries always leads to a smaller value of 77. If these binaries 
are not properly taken into account, the overestimation in the de- 
rived dynamical mass is given by Md yn /M c i = 9.75 77 J , where 
M c i is the true cluster mass. In our analysis we consider three 
types of clusters: 

- Particle-dominated clusters. The measured velocity disper- 
sion is dominated by the motion of the stars/binaries in the 
cluster potential. Clusters with a low binary fraction or a 
large stellar density (^ 10 7 stars pc -3 ) are good examples. 
In this case Eq. (1) applies and 77 « 9.75. 

- Intermediate clusters. Most realistic clusters are of this type. 
If the canonical value 77 = 9.75 is adopted, the derived value 
of Mdyn may mildly overestimate the true cluster mass, M c \. 

- Binary-dominated clusters. The measured velocity disper- 
sion is dominated by the orbital motion of the binaries, so 
that 77 <K 9.75. Eq. (1) may result in a significant over- 
estimation of Mdy n if the presence of binaries is ignored. 
Binary-dominated clusters generally have a low stellar den- 
sity (~ 0.1 stars pc -3 ) and a high binary fraction. 

The presence of binaries additionally affects the star cluster 
dynamics in the following way. Imagine a cluster consisting of 
N stars. Suppose we now add companions to each of these stars. 
The number of particles (i.e., singles and binaries) is still N, but 
the mass of the cluster has increased, resulting in larger centre- 
of-mass velocities than in the case of a single-star cluster. We 
will return to this issue in § 12. 41 



This paper is organised as follows. In §|2] we briefly de- 
scribe our method and assumptions. We discuss the effect of 
varying the star cluster properties, such as the size, stellar den- 
sity, mass, number of stars, virial ratio, stellar mass distribution, 
and the aperture size on the dynamical mass estimate in §[3] 
Subsequently, we study the effect of varying the binary popu- 
lation properties (binary fraction, mass ratio distribution, eccen- 
tricity distribution, and semi-major axis or period distribution) 
on the dynamical mass estimate in §|4] By varying each of these 
properties we study their respective contribution to the value of 77 
in Eq. (Q~|i. In §[5]we describe, from a practical point of view, un- 
der which conditions binarity can be ignored, and under which 
conditions ignoring the binaries results in a significant overesti- 
mation of the dynamical mass. Finally, we summarise and dis- 
cuss our results in §|6] 



2. Method and assumptions 

We study the effect of binaries on the dynamical cluster mass 
determination using simulated clusters. In our numerical simu- 
lations we can determine Md yn , cri os , and Rh m for each cluster, 
allowing us to derive the true 77 for clusters with different prop- 
erties, and to study the error that is introduced in Md yn when bi- 
narity is ignored (i.e., if the canonical value 77 = 9.75 is adopted). 

2.1. Model properties 

We simulate star clusters us ing the STARLAB package (see, e.g., 
iPortegies Zwart et al. I I2001I) . We draw N stars from a mass dis- 
tribution /m(M) (Mmin < M < M max ), such that the total 
mass of the cluster, after inclusion of the binary companions, 
is M c i = 10 4 M by default. We choose this value because it is 
a typical mass for open cluster-like objects, and because the ef- 
fects of binarity on the dynamical mass are pronounced for such 
clusters. Dynamical masses for such low-ma ss c lusters are de- 
rived i n, for example, Mandushev et al.l(ll99ll) and lde Grijs et al.l 
(2008) . Note, however, that most well-studied clusters have mea- 
sured dynamical masses of 10 5 ~ 7 M ; this is mainly because of 
(1) their brightness, and (2) their large velocity dispersion. We 
will discuss the full range of cluster masses in § !3.3l and §|5] We 
define the binary fraction as Fm = B/(S + B), where S and B 
are the number of single stars and binary systems in the clus- 
ter, respectively. The total number of "particles" is indicated by 
N - S + B (note that this is not the number of individual stars 
S + 2B). We ignore the presence of triple and higher-order sys- 
tems. 

A fraction Fm of the stars are assigned a companion star. 
The mass ratio, semi-major axis (or, alternatively, the orbital 
period) and eccentricity of each binary are drawn from the re- 
spective probability distributions f q (q), f a (a) and f e (e). We as- 
sume random orientation for the binary star orbits. Each particle 
(i.e., single star or binary) is given a certain position and veloc- 
ity according to the Plummer model (Plum merlll91 ll) . using the 
makeplummer routine in the STARLAB package. The cluster is 
scaled, such that it has a certain projected half-mass radius Rh m . 
For Plummer models, the corresponding intrinsi c half-mass ra- 
dius i s given by ^hm.intr — 

(2 2/3 _ i) i/2fl hm (e.g. JHeggie&Hul 
2003). The average mass density within the half-mass radius is 
given by (±M cl )/(f 7rfl 3 m ) = f^'M^ 3 . Similarly, the corre- 
sponding average stellar density is given by 

(p)hm = j?7i~ l (Mr) M c i R h ^ , (2) 
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where {Mj) is the average mass of a binary system. Each cluster 
is assumed to be in virial equilibrium, and no mass segregation 
is assumed to be present. 

2.2. Two canonical models 

In order to study the effect of each star cluster property sepa- 
rately, we perform our simulations with two different models, 
which we will refer to as (the simplified) model S and (the more 
realistic) model R. The default properties of the two models are 
listed in Tablefl] Throughout this paper we vary each property of 
the binary population individually, keeping the other properties 
constant, in order to study the effect of each binary parameter on 
77 individually. 

Model S is a simplified star cluster model, consisting of 
equal-mass stars. If the binary fraction in Model S is set to 
Fm = 0%, this model satisfies the assumptions of Eq. ([1]). We 
use this simplistic model, as changes of the cluster properties 
and the binary population have very pronounced effects on the 
derived cr los and rj, thus allowing us to quantify the relations pre- 
cisely. We assign a mass M = 1 M to each star, and adopt a 
binary fraction of Fm = 100%. We assign to each binary a mass 
ratio q = M2/M1 — 1, an eccentricity e — 0, and a semi-major 
axis a = 10 3 R (~ 4.65 au). 

Model R, on the other hand, is more realistic, and is a good 
approximation for real star clusters. Its properties are identical to 
those of model S, except for those mentioned below. Each star 
is assigned a m ass which is dr awn from the Kroupa initial mass 
function (IMF) (tKroupal2001l) , in the mass range 0.02 - 20 M , 
given by 



dN 

/Kroupa(M) = — oc 



M-° 3 for 
AT 13 for 
AT 23 for 



0.02 < M/M Q < 0.08 
0.08 <M/M <0.5 
0.5 <M/M <20 



(3) 

The Kroupa IMF is currently considered to be a good description 
of the stellar mass distribution. In our analysis we only draw 
stellar masses above the deuterium burning limit: M = 0.02 M , 
as lower-mass objects (such as planets) barely contribute to the 
dynamics of the cluster, the dynamical mass determination or the 
luminosity of the cluster. We adopt a maximum stellar mass of 
20 M Q . Stars more massive than 20 M Q are bright, but very rare. 
As we calculate the velocity dispersion directly from the line- 
of-sight velocities (i.e., we do not apply luminosity weighting), 
the results presented in the paper are practically independent of 
the choice for the the upper mass limit. The binary fraction for 
model R is 100%. We adopt a semi-major axis distribution of 
the form f a (a) oc aT x (10 R < a < 0.02 pc), a flat mass ratio 
distribution f q {q) = 1, and a thermal eccentricity distribution 
f e (e) = 2e. The latter choices are motivated in § 14.1114.21 and !4.31 
respectively. 



2.3. The velocity dispersion 

Direct measurements of the velocity dispersion of a star clus- 
ter can be done in three ways. Most commonly, cn os is deter- 
mined from (i) the width of spectral lines from observations 
integrated over a larg e part of the cluster (e.g., iBastian et al.l 
2006;|Moil et al.l[2007h . For nearby clusters the velocities of in- 
dividual stars can be measured, both (ii) radial velocities (e.g., 
Reims et al. 2006; Apai et al. 20 03) and (iii) proper motions 
(e.g., Ivan Leeuwen et all 120001 : IChen et all 120071: IStolte et all 



Table 1. The default properties of the two models used in our 
analysis: the simple model S (middle column) and the more re- 
alistic model R (right-hand column). In our analysis we vary the 
properties of each model in order to find the effect of this change 
on the value of 77. At the bottom of the table we list for each 
model the line-of-sight velocity dispersion <x los of the individ- 
ual stars, cr part of the centre-of-mass motion of the binaries, and 
cr 01 - b of solely the orbital motion of the binary components. Each 
value represents the width of the best-fitting Gaussian. 



Property 


Model S 


Model R 


Model 


Plummer 


Plummer 


Proj. half-mass radius 


R hm = 5 pc 


R hm = 5 pc 


Particles N = S + B 


N = 5 000 


N = 18 600 


Total mass 


M cl = 10 4 M Q 


M cl = 10 4 M 


Mass segregation 


No 


No 


Virial equilibrium 


Yes 


Yes 


Primary mass 


Mi = 1 Ms 


/Kro Up a(M 1 );0.02-20 M 


Binary fraction 


F M = 100% 


F M = 100% 


Mass ratio 


q=l 


/ ? (?) = 1;0< 9 <1 


Eccentricity 


e = 


/„(e) = 2e; < e < 1 


Orbital size 


a = 10 3 Ro 


/opi k (a); 10 R - 0.02 pc 


Orientation 


Random 


Random 


cr ios (kms -1 ) 


7.10 


1.20 


tr part (kms- 1 ) 


0.91 


0.91 


cr orb (kms^ 1 ) 


7.09 


0.14 



2007). Due to the nature of the observations, the velocity dis- 
persion obtained using techniques (i) and (ii) may be affected by 



the presence of binaries, while that obtained using technique (iii) 
is insensitive to binaries. Our paper thus applies to the spectral 
line and radial velocity studies, and not to proper motion studies. 

In our simulations we obtain radial velocities v r for each 
star in our simulated cluster. By fitting Gaussian profiles we de- 
rive the projected velocity dispersion <ti os from the projected ra- 
dial velocity distribution f Vr (v r ) of the stars in the cluster. We 
consider only radial velocities in the range (v r ) - 3 rms(vv) < 
v r < (v r ) + 3 rms(v r ), where (v r ) is the mean radial veloc- 
ity and rms(v r ) the corresponding root-mean-squared variation. 
Note that ~ 99.5% of the stars have line-of-sight velocities be- 
tween these limits. We use this cut-off in v r for the following 
reasons. The distribution over velocity is generally not Gaussian, 
and has broad wings (see Fig.Q]). The rare, extreme-velocity stars 
in these wings affect the determination of <xi os significantly. For 
example, the (dotted) fit in the bottom panels of Fig. Q] would 
be broader without the rejection of the 0.5% of the stars with 
extreme velocities. A similar cut-off is automatically imposed 
during spectral line analysis of unresolved star clusters, where 
the low signal-to-noise broad wings remain undetected, and the 
best-fitting Gaussian is essentially determined using the brighter, 
central part of the spectral line. 

The quantity cr\ os is affected by (i) the motion of the particles 
in the cluster potential, and (ii) the orbital motion of binary com- 
ponents about their centre-of-mass. We will refer to the orbital 
motion of the particles (i.e., the centres-of-mass) in a cluster as 
Cpan- We use cr or b to refer to the binary orbital motion, relative 
to their centres-of-mass. For the determination of the total clus- 
ter mass it is important to know whether one of these dominates. 
We refer to a cluster with cn os ~ cr part as particle-dominated, and 
to a cluster with cr\ os ~ cr orb as binary-dominated. For most clus- 
ters, however, both the particle motion and orbital motion are 
important; we refer to these as the intermediate case. Spitzer's 
equation is only applicable in the particle-dominated case, and 
results in an overestimation of Md yn in the intermediate and the 
binary-dominated cases. 
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Fig. 1. The contribution of the (centre-of-mass) particle motion 
and the (binary) orbital motion to the measured velocity dis- 
persion of a stellar population. The example shown here is for 
model R with = 5 pc and M c \ — 10 4 M and a binary frac- 
tion of 100%. Top panels: the line-of-sight binary orbital veloc- 
ity distribution / VoA (v rb) for model R, on a linear {left) and a 
logarithmic scale {right). Note that /„ orb (v or b) is independent of 
the cluster size, mass, and binary fraction. Middle panels: Same 
for the line-of-sight centre-of-mass particle velocity distribution 
/vpart( v part)> i- e -' on ly taking into account the motion of the centre- 
of-mass in the cluster. The particle velocity distribution is well 
fitted by a Gaussian distribution. Unlike the distribution in the 
top panels, this distribution depends strongly on the size and 
mass of the star cluster. Bottom panels: the measured line-of- 
sight velocity dispersion of the binary components in the simu- 
lated cluster. This distribution is a combination of / Vo[b (v OT b) an d 
/ Vpm (v part ). The best-fitting Gaussian is indicated with the dot- 
ted curve. The Gaussian function fits the central data well. The 
broad wings are not well fitted by the Gaussian; however, these 
are dominated by noise in real observations. 



For a binary population similar to that of model R (Table [TJ, 
the velocity dispersion resulting from orbital motion is <x 01 . b * 
0.14 kms" 1 Here we have included all binaries with equal 
weight, irrespective of their mass. This value is (by definition) 
independent of the binary fraction Fm- The weight given to cr 0lb 
relative to <x pal1 , when inferring cr los , however, depends on Fm- 
Note that the above value of cr 0lb is obtained using all stars in 
the cluster, although in reality the velocity dispersion is obtained 
for a very specific subset of stars. The latter selection effect may 



Table 2. When studying the effect of certain binary parameters, 
should we keep either the total cluster mass M c \, or the number 
of particles N = S + B constant? Below, we show the values 
of cn os and r\ for three different models. Model SI is identi- 
cal to model S, except that no binaries are present. Models S2 
and S3 have Fm = 100%. Model S2 has the same total mass as 
model SI, while model S3 has the same number of particles as 
model SI. Each single/primary star has a mass M = 1 M . Each 
binary has q = 1, e = and a = 10 3 R Q . Columns 2-4 list the 
number of particles {N - S +B), the binary fraction, and the total 
mass of each cluster. Columns 5 and 6 list the measured velocity 
dispersion cri os of the stars (i.e., singles and binary components) 
and of the centre-of-mass of the particles <x part . The correspond- 
ing values T] derived from these are listed in columns 7 and 8. 
As models SI and S3 have the same particle velocities cr part , we 
choose to keep M c \ (rather than AO fixed when comparing differ- 
ent clusters. Note that the values in column 8 are slightly larger 
than the canonical value of 9.75; see § l3.2l for details. 



1 


2 


3 


4 


5 


6 


7 


8 


# 


N 


Fm 


M a 


0"los 


°"part 


1 








% 


M G 


km s 


km s _I 






SI 


10000 





10000 


0.92 


0.92 


10.23 


10.23 


S2 


10000 


too 


20 000 


7.10 


1.29 


0.34 


10.24 


S3 


5 000 


too 


10 000 


7.06 


0.91 


0.17 


10.30 



change the effective value of <x orb significantly (see, e.g., § 14.21 
and §|43]i. 



2.4. Comparison issues 

In order to describe the effect of varying each binary param- 
eter on the dynamical mass derivation, a comparison between 
models "before" and "after" the modification is necessary, whilst 
keeping all other parameters the same. This requirement, how- 
ever, leads to some ambiguities. In particular, when adding bi- 
naries to a cluster, one has to make a decision whether to keep 
(i) the cluster mass M c \ or (ii) the number of particles N con- 
stant. By simply adding binary companions to a stellar pop- 
ulation, the average mass of each particle increases, resulting 
in a larger total cluster mass. The cluster mass is then given 
by M c i = N(Mi){l + (q)F M ), where (Mi) is the average pri- 
mary/single mass, Fm the binary fraction, M c \ the total cluster 
mass and {q) the average mass ratio (assuming a mass ratio dis- 
tribution that is independent of primary mass). Due to the larger 
mass, the particles move faster, such that ch os is larger. 

On the other hand, one could also decide to scale the num- 
ber of particles N in the cluster with binaries, such that its to- 
tal mass M c i is equal to a cluster without binaries. For a bi- 
nary population with a mass ratio distribution that is indepen- 
dent of primary mass, the number of particles N is given by 
N = M d (M 1 )- 1 {l + (q)F M r 1 . 

One thus has to make the choice to keep either N or M c \ 
constant. We use simulations to show the consequences of ei- 
ther choice. We perform our simulations with three models, 
S1-S3, for which the properties are listed in Table [2] Each 
model has a half-mass radius of 5 pc, and all stars have a 
mass of 1 M . Model SI is our reference model, which con- 
sists of N = 10 000 single, equal-mass stars, with a total mass of 
10 000 M . Models S2 and S3 include binaries. The binary frac- 
tion in these models is 100%, and each binary has a = 10 3 R Q , 
e = 0, and q = 0. Model S2 has the same number of particles 
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N as model SI, while model S3 has the same total mass M c \ as 
model S 1 . 

Model S2 has a larger velocity dispersion than model SI 
for two reasons, (i) the presence of binaries, and (ii) the in- 
creased total mass. In order to separate the effects of these two 
changes, we calculate cn os and 77 not only for the stars, but also 
for the (unmeasured) motion of the centres-of-mass of the bi- 
naries. Model S3 has a similar mass, but a smaller number of 
particles. The velocity of each particle, however, is very similar 
to that of model S 1 . 

Considering the motion of the centre-of-mass of each binary, 
it is better to study the effect of binarity using model S3 rather 
than model S2, i.e., to keep the total mass M c \ constant in the 
comparison, rather than the number of particles N. Unless stated 
otherwise, we therefore keep the cluster mass M c \ constant in 
each comparison, for the remainder of this paper. 

3. Dependence on cluster properties 

In this section we study the effect of varying the cluster prop- 
erties on the dynamical mass determination. We discuss varying 
the virial ratio in § 13.11 the mass distribution and aperture size in 
§ 13.21 the number of particles in § 13.31 and the half-mass radius 
in§E4] 



3.1. The virial ratio Q 

Spitzer's equation assumes that a cluster is in virial equilib- 
rium. Observations, however, suggest that many clusters form 
out of virial equilibrium (e.g.. iBastian & Goodwinl l2006). This 
often results in early dissolution into the field star population 
(in fant mortality) or significant mass loss (infant weight loss); 
see lde Gri js & Parmentieri d2007l) and references therein. 

The virial ratio Q = -Ejc/Ep of a star cluster is defined as 
the ratio between its kinetic energy Ek and potential energy Ep. 
Clusters with Q = 0.5 are in virial equilibrium, and those with 
Q < 0.5 and Q > 0.5 are contracting and expanding, respec- 
tively. Since Ek °c cr? , a star cluster has cr? oe Q. A more 
generalized version of Spitzer's equation, including the effect of 
virial equilibrium, is therefore 



M dyn = (20- 1 77 



los 



(4) 



which, if Q — 0.5, reduces to Eq. (HJ. If a cluster is assumed 
to be in virial equilibrium while in reality it is expanding, the 
dynamical ass overestimates the true mass by a factor of 2Q. 

The most important reason that many (if not all) clusters 
are formed out of virial equilibrium, is that the star-forming ef- 
ficiency is not 100%. After removal of the gas by the winds 
of the most massive stars, the gravitational potential of the 
cluster is r educed significantly, which results in cluster expan - 
sion (e.g.. iKroupa & BoiTvl 120021; IBastian & Goodwin! [2006). 
iGoodwin & Bastianl d2006l) define the effective star-forming ef- 
ficiency (eSFE) e as the star-forming efficiency that one would 
derive from the virial ratio under the assumption that the star- 
forming cloud was originally in virial equilibrium: Q = (2e) _1 . 
Under this assumption, the dynamical mass overestimation for a 
cluster of single stars is Md yn /M c i — 

3.2. The mass distribution and aperture size 

Spitzer derived Eq. ([l} for a cluster of equal-mass single stars, 
assuming that the line-of-sight velocity is measured precisely at 
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Fig. 2. The effect of the aperture size on cr los and 77, for model S 
(dashed curves) and model R (solid curves). Each model has a 
zero binary fraction, a half-mass radius Rh m = 5 pc, and a total 
mass M c i = 10 4 M . This figure shows the derived <xi os and 
77 for all stars within an aperture, as a function of the aperture 
size. The derived value 77 for each model depends only mildly on 
the mass distribution, while the aperture size is of much greater 
importance. Depending on the size of the aperture, the total mass 
of a star cluster may be overestimated by a factor of two, if this 
selection effect is not taken into account. 



the half-mass radius. This measurement is impractical to do; usu- 
ally an integration is performed over a large part of the cluster. 
In reality the outskirts of a cluster are dominated by the back- 
ground field star population, so that the analysis is restricted to 
the inner part of the cluster. These selective measurements in- 
troduce biases. The stellar density in the cluster centre is higher 
than that in the outskirts, resulting in a larger velocity disper- 
sion. The magnitude of this difference (and thus the difference 
in 77) depends on the density profile of the cluster. If this selec- 
tion effect is not taken into account, the derived Md yn depends 
on which part of the cluster is observed. 

Below, we study the dynamical mass overestimation as a 
function of position in two ways. First, we study the overestima- 
tion if the velocity dispersion is measured at a projected radius 
p from the cluster centre, as a function of p. Secondly, we study 
the dynamical mass overestimation for different apertures of in- 
tegration, as a function of the aperture size. The former can be 
calculated analytically for Plummer models, while the latter can 
be directly compared with observations. 

The line-of-sight velocity dispersion of the centres-of-mass 
at a certain projected distance p from the cluster centre is 



2 3nGM cl 
^ = 64 



1 + 



-,2 \ 



R lm> 



■1/2 



(5) 



for a Plummer model (e.g., Heg gie & Huill2003l) . At the half- 
mass radius, p = R^ m , this equation simplifies to 



0"part(tfhm) =2.116 



\ _1 



10 4 M, 



pc 



kms 



(6) 



Note that this is the projected velocity dispersion at projected 
radius p, rather than the integrated value within a radius p. The 
dynamical mass may be significantly overestimated if this se- 
lection effect is not taken into account. An expression for the 
dynamical mass overestimation as a function of p is obtained by 
substituting Eq. (0 into Eq. ([1]): 



M a 



V2 



1 + 



-1/2 



(7) 
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For velocity dispersions measured in the cluster centre this re- 
sults in a mass overestimation by ~ 40%. Measurements at the 
half-mass radius provide the correct Md yn , while measurements 
in the cluster outskirts result in an underestimation of the mass. 

In reality, the velocity dispersion of a star cluster is deter- 
mined from the measurements within an aperture. The size of 
the aperture is usually defined by the projected radius at which 
the projected stellar density becomes so low that the contribu- 
tion of background stars dominates the brightness and velocity 
dispersion. 

We illustrate the dependence of the derived dynamical mass 
on the aperture size by performing simulations of models S 
and R, each with a zero binary fraction, Fm = 0%. We use 
N = 10000 for model S and N = 28 000 for model R, so that 
each cluster has a mass of M c \ = 10 4 M . Fig. [2] shows <x los 
and T] as a function of the Lagrangian radius R f within which 
the measurement is performed. The (projected) nth-percentile 
Lagrangian radius Rf is defined as the radius which includes 
a fraction / of the total mass of a cluster. For example, the 
projected radius which contains 50% of the mass, Rso% (i.e., 
the half-mass radius Rhm), for example, is found by solving 
M(p)/M c \ = 50%. For a Plummer model, the mass M(p) within 
a projected radius p is given by 



M(p) = M c , 



1 + _hm 

P 2 



(8) 



(e.g.. lHeggie & Hull2003l) . The projected radius Rf which con- 
tains a fraction / of the cluster mass, is given by 



^hn 



(9) 



The choice of the aperture size may result in dynamical mass 
measurements differing by up to a factor of two. For example, if 
only the stars within a radius Z?20% = j^hm are observed, while 
this is not taken into account, the total mass is overestimated by 
~ 60%, for each model. The velocity dispersion of the model 
with the Kroupa mass distribution is slightly lower than that of 
the equal-mass model in the cluster centre, while the values are 
virtually the same for the cluster as a whole. The derived values 
of T] for the Kroupa model are 3% larger than those of the equal- 
mass model in the cluster centre, but practically equal if the en- 
tire cluster is taken into account. This indicates that Spitzer's 
equation is not very sensitive to the mass distribution. 



An open star cluster typically has a mass of 10 4 M©. 
Assuming a Kroupa mass distribution and a minimum mass of 
0.02 M , the number of particles (singles or binaries) in such 
a cluster is expected to be N ~ (1.4 - 2.8) x 10 4 . The up- 
per limit represents a cluster without binary systems, and the 
lower limit a cluster with Fm = 100% and a mass ratio of unity 
for all binaries. For young massive star clusters, with a total 
mass of order 10 6 M , the number of particles is expected to 
beN = (1.4- 2.8) x 10 6 . 

In order to study the relationship between N and the derived 
77 and Mdyn, we simulate models S and R, and evaluate the results 
for different values of N. Each model has a binary fraction of 
100% and a half-mass radius R^ m = 5 pc. The results are shown 
in Fig. [3] Due to the different mass and mass ratio distributions 
of each model, the average mass of a particle is different: 2 M in 
model S and 0.54 M in model R. Given the number of particles 
N, model S is therefore 3.69 times more massive than model R. 
This difference in total mass contributes to the larger velocity 
dispersion in model S, which is reflected in Fig. [3] 

According to Eq. (H), i] is proportional to Md yn o" lo ^ for a set 
of clusters with given half-mass radius Rh m . Since Md yn °c N, 
this can be rewritten as 77 oc Afov?. Eq. ([TJ applies for a particle- 
dominated cluster, irrespective of N. In a binary-dominated clus- 
ter, the measured velocity dispersion <x los is independent of N. 
For binary-dominated clusters we thus have 77 oc N. Ignoring the 
presence of binaries in a binary-dominated cluster therefore re- 
sults in a mass overestimation Md yn /M c i oc N~ l . Note that this 
relation is only applicable in the binary-dominated regime, and 
therefore cannot be extended to arbitrarily large N. 

These properties are clearly shown in Fig. [3] Note that, since 
we adopt a fixed half-mass radius of 5 pc, model S is only binary- 
dominated for N S5 10000. In this regime, model S shows the 
predicted constant cr\ os , and the corresponding 77 oc N. Model R, 
on the other hand, shows a stronger correlation between <x los and 
N; the overestimation of Md yn is also less severe. 

We finally describe the subtle difference between varying the 
number of particles N and the cluster mass M c \. These quantities 
are related as M c \ = N{Mj), where (Mr) is the average particle 
(single or binary) mass. When changing M c \ while keeping N 
fixed, one changes the average mass (Mr) of the binaries. The 
latter is equivalent to changing the mass distribution, the mass 
ratio distribution, or the binary fraction. We refer to § 13.21 § 14. 21 
and § !4.4l for a discussion on these issues, respectively. 



3.3. The number of particles N - S + B 

Under the assumptions made by Spitzer, Eq. ([TJ is independent 
of the number of particles. The only relation between Md yn and 
N has a statistical nature: if the number of particles is small, the 
statistical error on Md yn is large. If binaries are present in the star 
cluster, however, the derived dynamical mass decreases with in- 
creasing N, up to the point when cr los is dominated by the particle 
motion. This can be understood as follows. Imagine a group of 
binary systems, which have a given orbital velocity. This orbital 
velocity is independent of the number of stars in a cluster. The 
stellar density (which is proportional to N), however, affects the 
motion of the centre-of-mass of each binary system. The parti- 
cles move faster in a dense cluster, while the orbital motion of 
each star in a binary system remains the same. The contribu- 
tion of cr part to o"i os , relative to <x orb , thus becomes smaller; the 
dynamical mass overestimation is less severe for clusters with 
larger N. 



3.4. The half-mass radius R^m 

Given the number of particles N or the cluster mass M c \, the 
velocity at which the particles move is dominated by the size 
/?hm of the star cluster. Particles move fast in small (dense) clus- 
ters, and slowly in large (sparse) clusters. The orbital motion 
of binary components, however, is unaffected by the size of the 
cluster. In a large (sparse) cluster, we therefore expect the bi- 
nary orbital motion to dominate <x los , and in a small cluster the 
motion of the particles dominates. The dynamical mass overes- 
timation from Spitzer's equation is thus least severe for small 
(dense) clusters. 

We study the relation between /? n m, o"ios and 77 by simulat- 
ing clusters of different size. We adopt the properties of mod- 
els S and R, each with a binary fraction of 100% and a clus- 
ter mass of 10 4 M , and study clusters with sizes in the range 
0.01 pc < /?hm < 100 pc. The resulting trend of <x los and 77 with 
/?hm is shown in Fig. |4] Above each panel we indicate, {p\ m , the 
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Fig. 3. The effect of the number of particles N — S + B (horizon- 
tal axis) in the cluster, for model S (dashed curves) and model R 
(solid curves). Each model has a binary fraction of 100% and a 
half-mass radius R^ m = 5 pc. For a given number of particles 
N, model S is 3.69 times more massive than model R. In dense 
clusters (with large N) the motion of the binaries in the cluster 
potential dominates <x los , so that 77 is close to its "zero binary 
fraction" value of 9.75. In sparse clusters (with low N), <Ti os is 
dominated by the orbital motion of the binaries. In this case, the 
overestimation of Md yn is largest. 

p hm (stars pc" 3 ) for model R p hm (stars pc" 3 ) for model R 
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Fig. 4. The parameters o"i os and 77 as a function of the half-mass 
radius R^ m , for model S (dashed curves) and model R (solid 
curves). The dotted curve indicates the results for a cluster with- 
out binaries. Each model has Fm = 100% and M c \ = 10 4 M . 
Above both panels we indicate the average stellar density within 
the half-mass radius. As model S is binary-dominated, we find, 
to first order approximation, <t\ os m constant and 77 oc Z?r^. The 
effect of binaries is less pronounced for model R. 



average stellar density within the half-mass radius for model R, 
given by Eq. (0. 

The figure shows that, as expected, larger clusters have 
smaller values for <xi os and rj. The left-hand panel shows that 
o"i os for model S (dashed curve) is approximately constant for 
Rhm ^ 0.5 pc, indicating that such clusters are binary-dominated. 
For these clusters the orbital motion is independent of Rh m . In 
Spitzer's equation, 77 oc Z?rJ. which is indeed observed in the 
right-hand panel of Fig. [4] For clusters with Rh m S> 0.5 pc, on 
the other hand, cr\ os decreases with increasing Rh m , indicating 
that the centre-of-mass motion of the particles is not negligible 
as compared to the orbital motion; these clusters are of the inter- 
mediate case. 

The effect of binaries is less pronounced in model R (solid 
curves in Fig.@). The value of <ti os is affected by both the particle 
motion and the orbital motion, but mostly dominated by the for- 
mer. If cries were completely dominated by the particle motion 
one would expect cn os oc R^j 2 and 77 = constant. In Fig.|4]this is, 
to first order, the case, although the dynamical mass of clusters 



Fig. 5. The parameters cn os and 77 for models with different 
semi-major axis distributions f a (a), for model S (dashed curves) 
and model R (solid curves). Each model has a binary fraction 
Fm = 100%, a half-mass radius R^ m — 5 pc and a total mass 
M C ] = 10 4 M . The curves indicate the results for models where 
all binaries have identical semi-major axis f ao (a) = 5{a — ao), 
with ao indicated on the horizontal axis. In the binary-dominated 
case (ao < 10 au), we have, to first order approximation, 
77 oc a . The horizontal lines indicate the results for a distribu- 
tion /opik(a) °c cT x (Opik's law). The results for the log-normal 
period distribution /dm(P) (not shown) are practically indistin- 
guishable from those for Opik's law. 

with large R^ m is clearly overestimated, if the assumption of a 
zero binary fraction is made. For clusters with M c \ = 10 4 M Q 
and /?hm = 10 pc, for example, we find 77 w 5. This is half the 
canonical value of 77. Ignoring binaries in such a cluster may thus 
result in a dynamical mass overestimation by a factor of two. 

4. Dependence on the properties of the binary 
population 

In this section we discuss the influence of varying the different 
binary parameters on the dynamical mass determination. We dis- 
cuss variations of the semi-major axis (and period) distribution 
in § |4.1| the mass ratio distribution in § 14. 21 the eccentricity dis- 
tribution in § 14.31 and the binary fraction in § 14.41 

4. 1 . The semi-major axis distribution f a (d) 

The distribution of binary semi-major axes (or orbital periods), 
relative to the size of the cluster is one of the most important 
parameters affecting our interpretation of the observed velocity 
dispersion cr los . The orbital velocity of a binary system is pro- 
portional to gT 1 / 2 . 

In order to only extract the contribution of a to the mass 
derivation, we simulate clusters with varying semi-major axis 
distributions. We first study the results for clusters in which all 
binaries have an identical semi-major axis: f ao (a) = 6(a - ao). 
Each model has a binary fraction of 100%, a half-mass radius 
of 5 pc, and a total mass of 10 4 M . The dependence of cn os 
and 77 on the value of aq is shown in Fig. [5] for both models S 
and R. In a binary-dominated cluster cn os w cr orb , we therefore 
have cries oc cr 1 / 2 , and hence 77 oc a -1 , This is clearly shown 
for model S (dashed curves) in Fig. [5] The effect is less pro- 
nounced for model R, which is neither binary-dominated nor 
particle-dominated. 

In reality, stellar groupings contain binaries with a large 
range of orbital sizes, unlike the example described above. 
These orbital sizes can be quantified using a semi-major axis 
distribution f a (a), or, indirectly, an orbital period distribution 
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fp(P). The flat distribution in log a, commonly known as Opik's 
law, has been observed for a wide r a nge of stellar populations 
(e.g.. IQpik|[l924t Ivan Albadalll968t IVereshchagin et al.l ll988t 
iPoveda & Allenll2004tlPoveda et al.ll2007l) . and is equivalent to 

(10) 

with q m j n » 10 Rp an d a max ~ 0.02 pc (4500 au). 
Duq uennov & Mavorl d 1 99 lb studied binarity among solar-type 
stars in the solar neighbourhood and found a log-normal period 
distribution: 

f (logP-logP) 2 ) 

f DM (P) oc exp I - 2 Pmin <P< Anax , 

V log P ) 

(11) 

with P min * 4 days and P max * 0.3 Myr. They find logP = 4.8, 
0"io g p = 2.3, where P is in days. The latter distribution is often 
used as the standard reference for the orbital size distribution of 
a binary population. The size of the smallest orbits is determined 
by the radii of the stars or, more precisely, by the semi-major axis 
at which Roche lobe overflow occurs. The size of the largest or- 
bit is determined by properties of the surrounding stellar popula- 
tion, in particular the stellar density. In our models with the semi- 
major axis distribution f„{a) we adopt the limits a m i n = 10 R Q 
(* 0.05 au) and a max = 10 6 R Q (~ 0.02 pc). In the models 
with the period distribution /dm(P), we adopt P mm = 4 days and 
^max = 0.3 Myr. The results for Opik's law and the log-normal 
period distribution are indicated by the horizontal lines in Fig. [5] 
Both distributions show practically indistinguishable results. 

The upper limits a max and f max are in reality dependent 
on the environment, in particular on the stellar density (e.g ., 
iBahcall et alJl985tlClose et al.lfT990tlChaname & Gouid(l2 004). 
Given the half-mass radius Rh m and number of particles N in 
a cluster, there is a maximum semi-major axis at which a bi- 
nary system is marginally stable. Very wide binary systems are 
quickly ionised, as their binding energy is too weak to keep the 
components together. Below, we derive a simple expression for 
the maximum semi-major axis a max , relating it to N and Rh m . 

Consider a star cluster with N particles and a half-mass ra- 
dius /?hm- The particles within the half-mass radius occupy a 
total volume 1 7r/?jJ . The average volume available to one parti- 
cle equals §7ri?jLj> , which is equivalent to a sphere of radius 
(2Ny 1 ^ 3 Rhm- The maximum semi-major axis is thus expected 
to be of order a max w (2A0 -1 ^ 3 ^hm- Note that this is a conserva- 
tive estimate as the density in the cluster centre is significantly 
higher than at the half-mass radius; the maximum semi-major 
axis is likely smaller than a max . The semi-major axis distribution 
is further truncated for older clusters due to dynamical evolution. 

Our simulations show that the effect of our a max = 0.02 pc 
assumption is relatively small. We have compared the difference 
between clusters with a max = 0.02 pc, and those with the more 
realistic a max as (2N) ^ Rh m for a cluster with half-mass radius 
Rhm and N particles. The difference in the mass-overestimation 
is A(Md yn /M c i) * 4 - 6% for realistic values of a max . 

4.2. The mass ratio distribution f q (q) 

In this section we study the relation between the mass ratio dis- 
tribution f q {q) and the systematic error in Md yn caused by ne- 
glecting the binaries. For each model we first adopt a mass ratio 
distribution f qo (q) - 6{q — qo), and evaluate rj for < q < 1; 
see the top panels in Fig. [6] The figure indicates that o"i os in- 
creases with increasing go for model S. This effect is less promi- 
nent for model R. The value of rj decreases with increasing qo 
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Fig. 6. The values cr\ os and rj for model S (dashed curves) and R 
(solid curves), as a function of the mass ratio distribution f g (q). 
Each model has F M = 100%, R hm = 5 pc, and M d = 10 4 M Q . 
The curves indicate the results for models in which all binaries 
have a fixed mass ratio qo, which is indicated along the hori- 
zontal axis. The top panels show the results for all stars in the 
cluster. The bottom panels show the results for the stars more 
massive than 1 M Q (i.e., mostly primary stars). The horizontal 
lines in each panel show the results for the flat mass ratio distri- 
bution f q (q) = 1. For model R with a flat mass ratio distribution, 
the dynamical mass overestimation is Md yn /M c \ w 1 .7 if all stars 
are included in the fit, and Md yn /M c i « 3 if only the stars with 
M > 1 M are included. 



for both models, which is most pronounced for small mass ra- 
tios. For very small mass ratios rj reaches 9.75, the value for a 
cluster without binaries. The mass overestimation for model R is 
approximately constant for qo k 0.3. 

In reality, a cluster contains binary systems covering the full 
range < q < 1. We therefore also perform our simulations 
using a continuous mass ratio distribution f q {q). The mass ratio 
distributi on has been studied for vari ous stellar populations. For 
example, iKouwenhoven et ail (120071) find a distribution f q (q) oc 
q-0A f or i n termediate-mass sta r s in th e nearby OB association 
Sco OB2. iKobulnickv & Frverl (|2007l) find very similar results 
for Cyg OB2. Others find a large numb er of q 1 ("twin") bi- 
naries, mostly among massive stars ( e .g. , iPinsonneault & S tanek 
2006; Lucy 2006; Soderhielm 2007). In this paper we study the 
results for a flat mass ratio distribution fn- d t{q) = 1. This dis- 
tribution is solely chosen as an example, in order to show the 
(typical) effect of having a mass ratio distribution rather than a 
single mass ratio. The results for this mass ratio distribution are 
indicated by the horizontal lines in Fig. [6] For both models S 
and R, the mass overestimation is similar to that of a model with 
qo ~ 0.4; the distribution f q (q) = 1 can thus be described with 
an effective mass ratio g e ff ~ 0.4. 

In an unequal-mass binary system the most massive star or- 
bits with the smallest velocity. In general, the velocity dispersion 
measured from spectral lines reflects the properties of a specific 
subset of the cluster members only. These are often bright stars 
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Fig. 7. The results for different eccentricity distributions, for 
model S (dashed curves) and model R (solid curves). Each model 
has M cl = 10 4 M e , R hm = 5 pc and F M = 100%. The curves 
indicate the results for models with an eccentricity distribution 
fe (e) = 5{e - en), where en is indicated along the horizontal 
axis. The horizontal lines indicate the results for the flat distri- 
bution ,/fl a t(<?) = 1 and the thermal distribution /2 e (e) = 2e. Top: 
results for all stars. Binaries with a large eccentricity spend most 
of their time near apastron, where their velocity is relatively low. 
The dynamical mass derivation for clusters with a large average 
eccentricity is therefore less affected by the presence of binaries. 
For realistic models, however, the effect of varying the eccen- 
tricity on rj is small as compared to changes in the other binary 
parameters, such as the binary fraction and the semi-major axis 
distribution. Bottom: results for stars with masses greater than 
1 M Q (i.e., mostly the primaries). With additional selection ef- 
fect, the dynamical mass overestimation increases significantly: 

(Md yn /M c l)M>l M ~ 1-8 (Mdyn/Mcl)aU stars- 



with narrow lines. As the radial velocity amplitude of a primary 
star increases with increasing mass ratio, the value of cn os in- 
creases as well. The bottom panels in Fig. [6] illustrate this ef- 
fect. In this example only the stars more massive than 1 M 
(i.e., mostly primary stars) are included. We have selected a 
subset of more massive stars. The orbital velocity of binaries 

— 1/2 

follows v rb 00 M T . On average, the contribution of cr or b to 
0"i os is therefore larger. We are thus measuring a smaller 77 and 
derive a larger dynamical mass overestimation by looking at 
the more massive stars. Not taking this selection effect into ac- 
count may result in a further overestimation of the dynamical 
mass of the star cluster. In this example we find that by select- 
ing the stars with M > 1 M , the mass is overestimated by a 
factor of ~ 1.8 more than if we would have selected all stars: 

stars- 



4.3. The eccentricity distribution f e (e) 

A binary system with eccentricity e > spends most of its time 
near apastron, where the orbital velocity of both components is 
relatively small. The probability of the binary near periastron 



(i.e., of finding a large velocity) is small. For this reason, the 
derived value for the dynamical mass decreases with increasing 
average orbital eccentricity. 

In order to study the effect of varying the eccentricity on 
the derived dynamical mass, we perform simulations of mod- 
els S and R. We evaluate the results for three eccentricity 
distributions: the single-value distribution f eo {e) = 6(e - en), 
the thermal distribution = 2e, and the flat distribution 
/flat(e) = 1. The thermal eccentricity di stribution is expected 
from energy equipartition (Heggie1 [l975l) . This distribution is 
fre quently adopted in dynamica l models, for reasons of simplic- 
itv. lDuquennov & Mavorld 19911) show that observations of solar- 
type stars in the solar neighbourhood are consistent with fi e (e). 
For populations with a thermal eccentricity distribution, most bi- 
naries are in highly eccentric orbits (50% have e > 0.7). Binarity 
among these objects is particularly difficult to detect spectro- 
scopically, as these binaries spend most of their time near apas- 
tron. Although most observations are consistent with the thermal 
distribution, they are often equally consistent with the flat ec- 
centricity distribution /fl at (<?) = 1, which is the third distribution 
considered in this paper. 

The top panels in Fig. |7]show the results for models S and R. 
The curves indicate the results for f eo (e), where the eccentricity 
en is indicated along the horizontal axis. The increasing 77 for in- 
creasing average eccentricity is clearly visible. For any realistic 
star cluster (similar to model R), however, the contribution to the 
systematic error in rj due to variations in f e (e) is small as com- 
pared to that of the other binary parameters, such as f a (a) and 
Fm- For these clusters, those with e = are associated with a 
mass overestimation of ~ 80%, while those with e = 0.95 are 
associated with a mass overestimation of ~ 50%. The horizontal 
lines indicate the results for /2 f (e) and /fl a t(<?). For realistic mod- 
els (e.g., model R), the results for /2 f (e) and ff\ at (e) are very sim- 
ilar to those of the clusters with / e „(e): the corresponding mass 
overestimation is Mi ya /M c \ » 65 - 70%. 

The bottom panels of Fig.|7]show results for the same simu- 
lations, but now only for the subset of stars which are more mas- 
sive than 1 M (cf. Fig.[6]and § I4.21 >. The results for the distribu- 
tions ft e (e) and /fl at (<?) are again very similar to those of clusters 
with f ea (e). For each model, the mass overestimation is now 1 .8 
times larger: (M dyn /M c i) M>1 m ~ 1.8 (M dyn /M c i)aii stars ; cf. §|fk2] 
By selecting a subset of massive stars, the derived mass is thus 
significantly further overestimated. The reason that the mass se- 
lection effect results in a larger dynamical mass overestimation 
is that for the orbital velocities of binary stars v or b °c M T 1/2 
that the measured cn os is larger. 



so 



4.4. The binary fraction F M 

In this section we discuss the relation between the dynam- 
ical mass overestimation Md yn /M c \ and the binary fraction 
Fm- Observations and simulations have indicated that the vast 
majority of stars, possibly even all stars , are formed in bi- 
nary o r multiple systems (e.g., Mathieu 1994; Maso n~etail 
1 19981: iGoodwin & Kroupal I2005L iKobulnicky & Fryerl 12001 
iKouwenhoven et al. 2005, 2007, and references therein). The re- 
sults presented in the sections above are based on clusters with a 
binary fraction of 100%. These results are therefore reasonably 
accurate for young star clusters. The binary popula t ion in clus- 
ters is known to decrease with time: Sollima et al. (2007j) find 
Fm = 10 - 50% for 13 low-density gl obular clusters. The b inary 
fraction in dense cores is even lower: ICool & Boltonl (120021) . for 
example, find Fm < 7% for the globular cluster NGC 6397. 
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Fig. 8. The effect of the binary fraction Fm in the cluster, for 
model S (dashed curves) and model R (solid curves). The fig- 
ure shows that, as expected, the overestimation of Md yn increases 
with increasing binary fraction. Model S shows a strong increase 
in o"i os at Fm ~ 50%, while the effect for model R is approxi- 
mately linear. 



Clearly, if the binary fraction is negligible (Fm ~ 0%), the 
star cluster can be considered as a cluster of single stars. In this 
case the value of 77 for a single-star cluster can be used, and 
Eq. ([T) gives a good approximation to the true mass. The pres- 
ence of binaries gradually becomes more important if the binary 
fraction increases. When comparing two models with a different 
binary fraction, one has to make the choice between either keep- 
ing the total number of particles N = S + B constant, or keeping 
the total mass M^ yn of the cluster constant. For reasons described 



in § !2.4l we choose to do the latter. 

The results of varying the binary fraction are shown in Fig. [8] 
for model S (dashed curves) and model R (solid curves). For both 
models, the presence of binaries is relatively unimportant if the 
binary fraction is below ~ 20%. For larger binary fractions, bi- 
naries rapidly become more important. For model S, going from 
Fm = 40% to 60%, the binaries suddenly become important for 
cr los . The reason for this is that, for model S, the velocity dis- 
persion of the binary orbital motions is much larger than that of 
the particle motion. For a binary fraction 0% < Fm < 100%, 
the measured velocity distribution has two components: a broad 
component for the binary systems and a small component for 
the single stars. When fitting a Gaussian distribution to the (two- 
component) velocity distribution, the single stars dominate the 
fit for small Fm, while the binary systems dominate for large 
Fm- For non-linear least-square fitting of a Gaussian function, 
this results in a rapid increase of the best-fitting <x los when the 
binary fraction increases from 40% to 60%. For model R, on the 
other hand, the transition is much smoother, as the particle and 
orbital motions are comparable, and results in an approximately 
linear relation between Fm and 77. For model R in our exam- 
ple, the relation between the mass overestimation and the binary 
fraction can be approximated by Md yn /M c i as (1 - 0.36Fm) _1 . 
This is, however, not a general result; models with other proper- 
ties should be studied individually, as Md yn /M c i increases with 
decreasing M c \, decreasing N and increasing R^ m . 

5. When can binaries be ignored? 

In this section we provide general criteria for when the pres- 
ence of binaries can or cannot be ignored when determining the 
dynamical mass of a star cluster. We study clusters with differ- 
ent structural parameters (M c \, N, Rhm) and different binary frac- 
tions Fm- We keep all other properties fixed to the values for 
model R (see Table [TJ. The dependence of cr part on the clus- 
ter mass M c \ and the half-mass radius R^m is shown in Fig. |9j 



Particle dominated 
■n ~ 9.75 




o part =1 km s" 1 



Binary dominated 
r| » 9.75 



0.01 



0.10 1.00 10.00 100.00 
Half-mass radius R hm (pc) 



Fig. 9. When can the presence of binaries be ignored? The 
grey-shades in the figure indicate whether the derived dynam- 
ical mass is reliable (light) or overestimated (dark). The solid 
lines indicate models with a centre-of-mass velocity dispersion 
<T part = 1 kms -1 (bottom) and 10 kms -1 (top). In the bottom- 
right and central region, the dynamical mass of the cluster is 
overestimated by a factor > 2 and a factor of 1 .05 - 2, respec- 
tively. In the top-left region the dynamical mass overestimation 
is less than 5%. Each model has a binary population identical 
to that of model R, with a binary fraction of 100%. For models 
with a lower binary fraction, the lines move towards the bottom 
right. 



where we have adopted a binary fraction of 100%. The value of 
0" part as a function of M c \ and R m -n can be derived from Eq. (Q}: 
Cpart V-^ci/^hm- From dark grey to light grey, the colours rep- 
resent binary-dominated, intermediate, and particle-dominated 
clusters. The regions are separated by iso-dispersion contours 
of cr part = 1 kms -1 and 10 kms -1 , respectively. Clusters with 
these dispersions have an average density of ~ 50 stars pc~ 3 and 
~ 5 x 10 7 stars pc -3 within their half-mass radius, respectively. 
For a lower binary fraction, the lines move towards the bottom 
right of the plot. 

Fig.[10]shows from an observational point of view the results 
for clusters with a binary population such as for model R. These 
are valid for any Plummer-like model, irrespective of mass M c \ 
and radius R^ m . The horizontal axis in each panel represents the 
measured line-of-sight velocity dispersion <x los . The five curves 
represent models with a binary fraction of 0%, 25%, 50%, 75% 
and 100%. The left-hand panel shows the contribution of the par- 
ticle motion to cr los . Clearly, the difference cr los - cr part increases 
with increasing binary fraction. Although this difference appears 
small, in the range of 0.1-0.5 kms -1 , its effect on the derivation 
of the dynamical mass can be large. The middle panel shows the 
derived value of rj as a function of cn os , for the different binary 
fractions. For large <x los we find 77 * 9.75 (dotted curve), indicat- 
ing that Spitzer's equation is a good approximation. The right- 
hand panel of Fig. [10]shows the dynamical mass overestimation 
Mdyn/Md. The mass overestimation increases with increasing bi- 
nary fraction and decreasing <xi os . The panel shows that if one 
measures cn os =6-8 kms -1 , the dynamical mass overestima- 
tion due to binarity is expected to be 10 - 30%. The right-hand 
panel of Fig. [10] additionally shows the expected deviation of 
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Fig. 10. When can binarity be ignored? This figure shows from an observational point of view how the measured velocity dispersion 
o"i os should be interpreted (adopting a binary population as in model R). The solid curves in each panel represent models with a 
binary fraction of 25%, 50%, 75% and 100%, respectively. The dotted curves indicate the results for a cluster with zero binary 
fraction, and are also representative of the particle motion in each cluster (irrespective of the binary fraction). Left: the difference 
between the measured line-of-sight velocity dispersion cn os and the intrinsic centre-of-mass velocity dispersion cr part . Middle: the 
true value of r\ as derived from the simulations. Right: the dynamical mass overestimation under the assumption that no binaries 
are present, i.e., using tj = 9.75. On the right axis of this panel we additionally indicate the shift in log(Ly/Md yn ) introduced by 
this systematic error. This figure is applicable to all clusters with a Plummer density distribution and a binary population such as in 
model R, irrespective of their size, mass, or central density. 



the datapoints in the log(Ly/Mdv n ) vs. age dia gram (see, e.g., 
iBastian et al .120061: lde~Griis & Parmentierl2007l) . due to binarity. 
Here, we have used the simple calculation Alog(Ly/Md yn ) = 
log(Ly/M d ) - log(L v /M dy „) = log(M dyn /M d ). Clusters with 
cr los > 3 kms -1 have Alog( 'L v /M dyn ) = 0.05 - 0.10. For clus- 
ters with o"i os < 3 kms" 1 , the inconsistency can be significantly 
larger. 

Most clusters in the literature for which the dynamical mass 
is obtained using Eq. (Q~|i are rather massive, often of order 
M c i = 10 5 " 7 M . This is mainly because of two reasons: (i) they 
are bright, and thus easy to detect, and (ii) they have a large ve- 
locity dispersion which is easily measured using high-resolution 
spectroscopy. As the observed clusters are generally massive and 
have a large cn os , their dynamical mass overestimation due to bi- 
narity is expected to be mild, of order 5%. 



6. Discussion and conclusions 

The total mass of a distant star cluster is often derived from 
the virial theorem, using line-of-sight velocity dispersion mea- 
surements and half-light radii. This dynamical mass, Md yn , is 
given by Spitzer's equation, Eq. (Q]i, under the assumption that 
no binary or multiple systems are present. This assumption is 
frequently made in the analysis of star clusters, for reasons of 
simplicity, although it is known that most stars are part of a bi- 
nary or multiple system. Ignoring this fact may lead to a signif- 
icant overestimation of the cluster mass. In this paper we have 
studied the validity of this assumption, and how this affects the 
dynamical mass determination. 

The measured line-of-sight velocity dispersion cr\ os of an un- 
resolved star cluster can be derived from spectral line analysis. 
The value cr los represents that of the individual stars. The veloc- 
ity of these stars is not only determined by the binary centre-of- 
mass motion cr part in the cluster potential, but also by the orbital 
motion cr or b of the binary stars. If the assumption is made that 
no binaries are present, cr\ os overestimates the motion of the bi- 
naries in the cluster potential. Application of Spitzer's equation 



without taking this effect into account then results in an overes- 
timation of Md yn . 

Whether or not binaries are important depends on (i) the 
star cluster properties: (ii) the properties of the binary popu- 
lation; and (iii) observational selection effects. In our analysis 
we therefore distinguish three cases: (a) the particle-dominated 
case, where binaries can be neglected and Spitzer's equation is 
approximately valid; (b) the intermediate case, for which Md yn 
is overestimated by no more than a factor of two, and (c) the 
binary-dominated case, where the orbital motion of the binaries 
dominates cn os , and Md yn is significantly overestimated. Particle- 
dominated clusters have cri os > 10 kms" 1 , intermediate clusters 
have 1 ^ <r los ^ 10 kms" 1 , and binary-dominated clusters have 
<x los ^ 1 kms" 1 . Depending on the cluster properties, binarity 
can introduce a shift of A log(Ly /Mi y n) = 0.02-0.5 in the mass- 
to-light ratio vs. age diagram that is often used to study star clus- 
ter evolution. The exact values depend on the properties of the 
star cluster and its binary population. 

The dynamical mass of a star cluster is overestimated if bi- 
nary systems are present but not properly taken into account in 
the analysis. Whether or not this overestimation is negligible de- 
pends on: 

(a) Star cluster properties. For a star cluster with a non-zero 
binary fraction, the structural properties of the cluster are impor- 
tant for the derivation of Md yn - Star clusters with a larger stellar 
density are least affected by the presence of binaries. These are 
clusters with a small half-mass radius Rh m , a large number of 
particles N = S + B, or a large total mass M c \. Spitzer's equa- 
tion assumes that the star cluster is in virial equilibrium, i.e., the 
cluster has a virial ratio Q = 0.5. If this is incorrect, the cluster 
mass is overestimated by a factor 2Q. 

(b) Binary population properties. Clearly, the dynamical 
mass of star clusters with a larger binary fraction Fm is more af- 
fected by the presence of binaries. For a binary-dominated clus- 
ter, binarity suddenly becomes important for Fm ^ 50%, while 
for the intermediate case the transition is gradual. Whether a 
cluster is particle-dominated, of the intermediate case, or binary- 
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dominated, is mainly determined by the ratio a e ff/^hm, where a e ff 
is the typical binary semi-major axis. In the binary-dominated 
case we have 77 oc Na^R^. Models with a smaller average 
eccentricity and a larger average mass ratio tend to be associ- 
ated with a larger dynamical mass overestimation. Variations due 
to the uncertainty in f e {e) and f q {q) are rather small, however: 
A(M dyn /M cl ) ^ 5%. 

(c) Selection effects. Observational selection effects are im- 
portant for the dynamical mass determination. Firstly, the dy- 
namical mass overestimation is more severe if the analysis is 
more concentrated in the central region of the cluster; this 
may result in a dynamical mass overestimation of up to 40%. 
Secondly, the velocity dispersion is often determined from spec- 
tral lines. The flux of these spectral lines is often dominated by a 
very specific group of stars, e.g., stars of a specific spectral type. 
If this specific subset of stars is not representative for the cluster 
was a whole (e.g., more massive stars are ofte n found in close, 
equal-mass binaries; IZinnecker & Yorkell2007l) . M dyn will be bi- 
ased. In § |4.2| and |4.3| we illustrated that our results differ if we 
study the velocity dispersion for a different stellar mass range. 
We find that by selecting stars with M > 1 M in our example, 
the mass overestimation increased by 80% with respect to the 
unbiased sample: (M dyn /M d ) M >i m q ~ 1-8 (M dyn /M c i) a n stars . In a 
subsequent paper, we will thoroughly study how such selection 
effects affect the derivation of M dyn . 

We have adopted several other assumptions in our analy- 
sis, which require further investigation. For example, we have 
assumed that the properties of the binary population, Fm, Mi, 
q, e, a (or P), the orientation, and orbital phase, are uncorre- 
lated. We made this assumption because the purpose of our pa- 
per was to illustrate the general effect of binaries on the dy- 
namical mass determination. Observations, however, have in- 
dicated that correlations between several of these parameters 
may be present. The binary fraction, for examp le, is known to 
decrease with decreasing prim ary mass (e.g., iFischer &"M arcv 
Il992t ISterzik & Durisen|l2004l) . and the eccentricity e decreases 
with decreasing orbital period P due to circularization (e.g., 
lHalbwachs et alj|2005l) . In the spectral analysis of a star cluster, 
one is generally sensitive to only a specific stellar mass range. 
These often bright or evolved stars may exhibit peculiar proper- 
ties that are not representative of the cluster as a whole. When 
interpreting specific observations, one has to specifically charac- 
terise the binary population of the sampled subset of stars, before 
deriving the dynamical mass of the cluster. 

Furthermore, we have not included triple or higher-order 
systems in our simulations, as the properties of these systems 
are very poorly const r ained by observations. The s tudies of 
iTokovinin & Smekhovl d2002) and lCorreia et alj d200fj|) find that 
20 - 30% of the wide visual bin aries have a spectroscopic sub- 
system (see also iHu et alj|2008h . However, our results are un- 
likely to be affected by this assumption, as stable multiple sys- 
tems are necessarily hierarchical. The vast majority of triple sys- 
tems, for example, consists of a primary star, a close spectro- 
scopic companion with semi-major axis a\, and a wide visual 
companion with semi-major axis 02- In or der for the triple sys- 
tem to be stable, we must have a 1 <k 02 (see lMardling & Aarsethl 
1200 1L for a detailed study). For the orbital velocities v\ and V2 of 
the inner and outer stars we have, to first order approximation, 
v\/v2 °c Vfl2/fli s> 1. Non-inclusion of the outer component 
would result in a mild additional overestimation of cr part , and 
thus of M dyn 

There are several other causes which may lead to a dynami- 
cal mass overestimation. Due to the stochastic and fractal nature 
of star formation, young star clusters are often found in groups. 



Confusion may lead to a measured velocity dispersion that is 
dominated by the systemic velocity difference between two star 
clusters, rather than the particle or binary motion of the systems 
in each cluster (e.g jMoll et alj|2007l) . 

Another bias in the dynamical mass may be caused by mass 
segregation. The mas s segregation can be prim ordial or dynami- 
cal (e. gjHunter et alJll997HBonnell & Davieslll998 l; IChen et all 
120071; iMcMillan et alj|2007l) . iBoilv et all d2005l) . for example, 
find that for clusters with a projected density of ~ 10 4 M Q pc~ 2 
or more, the assumption of r\ — 9.75 results in a significant dy- 
namical mass und erestimation. The latter result is confirmed by 
iFleck et ail d2006l) : as heavy stars quickly sink to the centre of 
the cluster potential, and as these dominate the spectral lines that 
are used to obtain cri os , the dynamical mass may be und eresti- 
mated by a factor of two. We refer to IFleck et ail d2006l) for a 
detailed discussion of the effect of mass segregation on 77 and 
M dyn . 

Although we have made several assumptions in this paper, 
our main conclusion is robust: the presence of binary systems in 
a stellar population often results in a significant overestimation 
of the dynamical mass, if one applies Spitzer's equation without 
properly taking into account multiplicity. 
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